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We study a one-dimensional class of reaction-diffusion models on a 10— parameters manifold. The 
equations of motion of the correlation functions close on this manifold. We compute exactly the 
long-time behaviour of the density and correlation functions for quenched disordered systems. The 
quenched disorder consists of disconnected domains of reaction. We first consider the case where 
the disorder comprizes a superposition, with different probabilistic weights, of finite segments, with 
periodic boundary conditions. We then pass to the case of finite segments with open boundary 
conditions: we solve the ordered dynamics on a open lattice with help of the Dynamical Matrix 
QO ' Ansatz (DMA) and investigate further its disordered version. 
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I. INTRODUCTION 



In this work we investigate in detail some exact properties, in ordered and disordered media, of a wide class of 
one-dimensional single-species reaction-diffusion models n which play a central role in the description of interacting 
many-particle systems in physics, chemistry, biology, etcEl . We consider systems on the 10 — parameter manifold on 
which {in arbitrary dimensions) the equations of motion of correlation functions close, and the dynamics is soluble. 

So far, most theoretical studies on Reaction-Diffusion models consider infinite, ordered media (see e.gcTta and 
references therein). We show that for times that exceed the typical diffusion time needed to span the size of the 
sample, i.e. in the regime where L 2 /t <C 1, (where L denotes the size of the lattice chain and t the time, see below), 
finite-size effects dramatically affect the dynamics. A (quenched-) disordered system consisting of a random collection 
of reaction-domains (each segment, of length C weighted with a probablity lo(C)) of varying sizes, when u>(C) is a 
decreasing exponential, the density exhibits a stretched exponential decay in the diffusion-limited pair annihilation 
C^I ■ and diffusion-limited coagulation mpdelsEl. This exponential decay also applies to the autocorrelation function of 
symmetric exclusion process (SEP)B. The picture is realistic for media, where disorder fractures the medium into 
disconnected reaction zones. The stretched exponential decay of the density of the photogenerated excitons in (MX) 
and (TMMC) chains is an experimentally established phenomenon!! For other distributions u)(C) and other reaction- 
diffusion systems the situation is more involved. However, in the regime where L 2 jt <C 1, the gwenc/ied-disorder always 
• affects the dynamics. Below, we study the density and the two-point non-instantaneous correlation functions for a 
t-H | class of reaction-diffusion systems in presence of quenched disorder. 

This paper is organized as follows: in the first part we study reaction-diffusion processes occuring on the soluble 
10— parametric manifold on which the equations of motion of the correlation functions are closed. In section II 
we determine this manifold as well as the equations of motions of the correlation functions. In section III. A we 
obtain the density and non-instantaneous two-point correlation functions for a finite and periodic chain. In section 
III.B, via the Dynamical Matrix Ansatz (DMA), we generalize previous works by Stinchcombe and Schutzcl, on the 
soluble 10— parametric manifold and obtain original results. The approach developped allows the exact computation of 
correlation functions on an open chain of arbitrary length, with injection and evacuation of particles at the boundaries. 
O ■ In the thermodynamic limit we compute the exact density and the non-instantaneous correlation functions. As a by- 
product we generalize the "Bethe Ansatz" equations previously obtained by Stinchcombe and Schiitz for the isotropic 
chain with symmetry-breaking boundary fields. In section IV, we introduce the disordered systems which we consider 
in the sequel. In section V, we take advantage of the results of section III. A to study the dynamics of quenched- 
disordered systems, where the quenched-disorder consist in a superposition of periodic segments with probabilistic 
weight oj(C). We consider the cases where to(C) is exponential and where it is an algebraic function. In section VI 
we use the original results obtained in section III.B to study the dynamics of quenched-disordered systems, where the 
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quenched-disorder consist in a probabilistic superposition of open segments with probabilistic weight u>(C). In section 
VII, we illustrate previous discussions by considering two physically motivated examples. 



II. THE FORMALISM 

It is known that models of stochastic hard-core particles are soluble on some manifold on wich the equations of 
motion of their correlation functions closdj. 

Consider an hypercubic lattice of dimension d with N sites (N — L d ), where L represents the linear dimension 
of the hypercube. On the lattice, local bimolecular reactions between single-species particle A, with hard-core, take 
place. Each site can be empty (denoted by the symbol 0) or occupied at most by a particle of type A denoted 
in the following by the index 1 . The reactions are specified by the transition rates V^L, where a, (3,^,5 = 0,1: 
V(a,/3)^(7,*),r2j:a + /9— ►t + J. 

Probability conservation implies T™^ = — J2( a P') ^"p anc ^ ^ap — ^' ^( a > P) ^ (T) ^) The state of the system 

is represented by the ket \P{t)) = J2{ n } where the sum runs over the 2 N ,N = L d configurations. At 

site i the local state is specified by the ket \m) = (1 0) T if the site i is empty, \rn) = (0 1) T if the site i is occupied 
by a particle of type A (1) . 

We define the left vacuum (x\ = J2{ n }({ n }\> which local representation is (\\ =(11)® (1 1). 

It is by now well established that a master equation can be rewritten formally as an imaginary time Schrodinger 
equation: Jr|-P(t)) = — H\P(t)), where H is the Stochastic Hamiltonian which governs the dynamics of the system. 
In sflHeral, it is neither hermitian nor normal. Its construction from the master equation is a standard procedure (see 
e.g.lilfi). The evolution operator H acts locally on two adjacent sites, H = J2 m Hm,m+i and because of .-probability 
conservation, we have (x\H = . The explicit form of the H considered here can be found, e.g., in Ref.cHa. 

Below we shall assume an initial state |_P(0)) and investigate the expectation value of an operator 0:(0)(t) = 
(x\Oe~ Ht \E(0)). For general single-species bimolecular reaction-diffusion systems, there are 12 independent 
parametersa. If one imposes to these parameters the 2 following constraints: 

pXOiT-ill /t-i00 , t-iOIA _ -plO i -pll ( r O0 , t-i01\ . pOl.-pll / r 00 , t-i10\ _ -pOl , -pll /-pOO , fi\ 
1 00 "+" 1 00 — \ L 11 + 1 11 ) — 1 01 "+" 1 01 ~ \ L 10 + 1 WJ ! 1 00 + 1 00 — \ L 11 + 1 11 ) — 1 10 + 1 10 ~ \ L 01 + 1 01J l W 

the equations of motion of the correlation functions close and the system is formally soluble (in arbitmaj dimensions) . 
It is useful to point out that these solvability constraints allow, through a similarity transformation!.™, to transform 
the equivalent quantum chain , into another ouantum chain sharing the same eigenspectra of an XX Z quantum chain 
with surface fields, of the type solved in Ref.L3. 

III. DENSITY AND NON-INSTANTANEOUS CORRELATION FUNCTIONS FOR ORDERED CASE, IN 

ONE-SPATIAL DIMENSION ON A FINITE LATTICE 

For the sequel, in order to compute the sampling average of quantities such as the density and the non-instantaneous 
correlation functions we need to know the expression of these quantities on a finite lattice. In the following subsection 
we briefly derive the latter expressions on a periodic finite lattice. In the subsection III.B we compute, as original 
results, the same expressions on an open lattice with help of Dynamical Matrix Ansatz. 

A. The periodic case 

For the class of systems which are described on a 10-parameter manifold^, according to the solubility constraints 
(Q), the equation of motion of the density, for a periodic chain of C < L sites, in one-spatial dimension, readsa: 

^W(t) = j t (x\n m e- Ht \P{0)) =A + B(n m )(t)+C(n m+ i)(t) + D(n m . 1 )(t) (2) 



1 This soluble 10— parametric manifold includes stochastic systems such as Glauber and voter models, the Symmetric Exclusion 
Process (SEP). 
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where we have defined: 

a — of 11 _i_ r!0 _i_ rOl. n — /of 11 _l rio , -pOi , -poo , -poo , -pOi , t-iio\ 

/1 = zi OO + i OO' i OOi - D = ~l 1 zi OO + i OO + i OO' i 10' i 01' i 10' i Olj 

(7 = rS; + r5i-(r55 + r5g); D = r?J+r^-(r;j + r^) ; e = ±Vcd-, m = ±/§ (3) 

In this work we focus on the case where E > (which also implies /i > 0). 

Solving the equations of motion obtained from we obtain in the ordered case (for integer m and m! and for 

•periodic segments of length £ < L): 

(n m (t)) c = ^£ m , Epei.B.z. (K'(0)>£ - (n m /(oo)) £ ) M m - m '(si5n J Er- m ' cos [p(m - m')] e 2 ^'™ 8 *', 

where (n m '(oo)) = 0(1 — <5a,o)(1 — fe+c+r>,o), with = — b+c+d » an< ^ the sum m ' runs over a ^ the integers 
belonging to the segments of length £ under consideration. The sum over p runs over the first Brillouin zone (I.B.Z.): 
P=^T, n = 0,...,£-l. 

With the same notations, we have for the non-instantaneous two-point correlation functions: (n m (t)ni(0))c = 
*TT Em' Epei.B.z. ((n m >(0)ni{0)) c - (n m > (°c)m{0)) c) fi"' {signE)"' cos [p{m - m')] e 2 \ E \ tco ^, where 
(n m '(oo)n.;(0)) = (n m '(oo)). 

It is useful for the sequel to obtain the long-time behaviour of these quantities in the regime where C 2 /\E\t 
L 2 /\E\t ~ 1, with C « L > 1 and \E\t > 1, for the density we have : 

(n m )(t) - (n m )(oo) = J2» m ~ m (M0)k - (V(oo)k) — == (4) 

, 2\/TV\E\t 
m' v i i 

The long-time behaviour (C 2 /\E\t ~ u = L 2 /\E\t ~ 1, with £ « L » 1 and \E\t 3> 1) of the two-point correlation 
functions reads 

(n x (t)n xo (0)) = E ^ V \ I — tItmT — — — ((%(0)n ao (0)) £ - {n y (^)n XB (0))c) (5) 



B. Density and non-instantaneous correlation function for the ordered case with open boundary conditions 

It is instructive to consider the physically motivated case of a disordered lattice with open boundary conditions. 
To investigate this situation, we begin with an ordered lattice of length L and open boundary conditions: we assume 
that at site 1 particles can be injected in the system with rate 5 and evacuated with rate 7. At site L particles can 
be injected in the system with rate 5 and evacuated with rate (3. We apply the Dynamical Matrix Ansatz (DMA) 
introduced.by Stinchcombe and Schutz to study the ordered Symmetric Exclusion Process (SEP) with open boundary 
conditional. n 

For two-states models, the Dynamical Matrix Ansatz assumes that the probability can be encoded byB \P(t)) = 

((W\ {Uj=i(M(t) +D(t)ar)\ ft)} |F))^, where Z L = ((W\C L \V)) is the normalization constant, C = E{t)+D(t), 

and o~~ denotes the usual Pauli matrix acting on site j. We have introduced the vectors ((W\ and |V)) on which 
act the (infinite and time-dependent) matrices E_(t) and D_{t) (not to be confused with the combination of rates (0)). 
To investigate the dynamics, we associate a spin-down to a particle and a spin-up to a vacancy. In so doing, the 
ferromagnetic ground-state | ft) corresponds to an empty lattice. 
In this formulation, correlation functions read: 

{n n (t)n n {t) . ..n jm (t)) = {{W^' 1 D(t)C^~ l D{t) . ..C L -^\V))^, j m > j m -x > ••• > k- 

The master equation, in its Hamiltonian formulation, reads + ff) \P(t)} = 0, where H — Y2j=i Hjtj+i +^1 
The first term in this expression corresponds to the bulk contribution, while the terms b\ and &l correspond to the 
boundary at site 1 (injection with rate a and evacuation with rate 7 ) and L (injection with rate (3 and evacuation 
with rate 5 ), respectively..—. 

Proceding as in RefBHO, on the 10— parametric manifold ([!]), we obtain the bulk contribution to the master 
equation: 

= ((7i + a A )C - (<7 3 + a 5 )D{t) + a 2 CD_(t)C_- 1 + aeC^D^C (6) 
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S{t) = (d! - a 4 )C + (ct 5 - <r 3 )D(t) + aaCDj^Cr 1 - (T 6 Cr l D{t)C (7) 

a 7 D(tf(t) = anCD(t) + a 9 D(t)C - T^C 2 + a 2 CD,(t)(r 1 D(t) + aeDfflCT 1 D(t)C, (8) 
Similarly, the boundary terms give rise to the following equations: 

0=((W\ {(2cr 4 - a)C- (2<t 5 -a- j)D(t) + 2a 6 CD,(t)CT 1 } = {(2oi - 8)C - (2<r 3 -0-S)D(t) + 2a 2 CD_{t)Cr 1 ^ \V)), 



where 



<n = rll + r£° > o ; a 2 = r£° + r£ - o x = c ■ a 3 = r°° + r°J + a x > o 
^ = rJJ + r°J > o ; a 5 ee r°° + rj° + * 4 > o ; <j 6 = r% + r}J - a 4 ee d 

oi = cr 3 + fT 5 - (r?5 + r?J + rj" + rJJ + rJJ + rJJ) ; cr 8 ee o-i + rJJ -rj} ; o- 9 = o- 4 + rJJ -rJJ, (9) 



with era = C, erg = D, (7i +ipj^-j= ^4 and (73 + (75 = — £J, where A, B,C, D are the quantities defined in (H). For the 
h&m, we have rj? = T% = _L_^ 



SEP model, considered in&ffiL3, we have Tq" = T^q = p f„_^ and thus with 172 = (73 = (75 = (76 = = T and 
(7i = (74 = as = <7g = 0, we recover previous results Refi 

To study the dynamical equations (^) with boundary conditions (|9j), it is useful to introduce the fol- 
lowing notation: 22 -(i) = Ql 22(*KZ _J , and so the correlation functions read (n,j 1 (t)rij 2 (t) . . . rij m (t)) = 
({W\D ji (t)D j2 (t)...Dj m (t)C L \V})^-, j rn > j m -i > ... > ji. For the study of the dynamics of the system 
we can then introduce the Fourier transform of D_{t)j{t) (for p ^ 0): T> p (t) = J^j e%V3 D-j(t)- Seeking a solution of the 
dynamical equation (||) (for p ^ 0) of the form V p (t) — e~ fpt V p (t — 0), we find e p = (03 + (75) — a 2 e~~ %p — <7 6 e lp . The 
equation (||) implies the relation: 2? P1 (t)2? P2 (i) = S P2tPl T> P2 (t)T> Pl (t), where S P2tPl is the amplitude of the scattering 
matrix which plays a central role in the theory of integrable models. Here, we have 

s P2, P i = - — - — i( V ,+n 2) i^r ' (pi>Pv r 0,0) (10) 

\(76 + a 2 e l ( pi+P2 > _ cr 7 eipi j 

For the SEP model we recover!'! the scattering amplitude of the isotropic Heisenberg chain 
With the definitions, 

5-2(74 - /3 + 2(q 1 -(7 3 ) _ 7 + 2((7 4 -g 5 ) of _ ?-2ai 

2a = ^^' 2/3 = 1 + — 2^ ;2t = i + — ^ — ' 2(5 = ^T' 

it follows from the boundary conditions that 



T(p) ee ((W\V P C L \V)) = ^"'^ e 2 *T(-p) = - Mf> S K 2 * L T(-p), (p ^ 0), 



(12) 



where f p (a,b) = 2(a + b) — 1 + e lp . These equations have as solution: T(p) = X)/>o a ' ( e ' pZ ~ eIf> ^ ^ 7 P (i"^ 

Consistency implies the following relation which quantizes the momenta p's, as in 



E;>o a ' (' 



quantum Many-Body problems: 



/p(a>7)/p(P>o) 



and more generally, for the m— points correlation functions, we have e 2lp ^ L 1 ' = B Pj Y\ 7 ^ =1 ^ k _ i -- ) ■ Equa- 

tions (|l3U.rand its to— points generalization are the analog of the Bethe Ansatz equations for quantum Many-Body 
problemsa'S. 

Correlation functions are then obtained as follows: 

(n n (t) . . . n jm (f)) - (n n (00) . . . n jm (00)) =/.../ II™ 1 [dpie-**-^*] T(pi, . . . ,p m ), 
j m > . . . > ji ; p 4 7^ ; 1 < i < m. Where T(pi, . . . ,p m ) ee ((^1^(0) . . ■'D Pm (0)\V))-^- encodes the initial state of 
the system. 

For finite systems the integration over the momenta, should be carried over p's satisfying the relations (|l3|). For 
infinite system (when L — > 00) these relations (Oh can always be fulfilled (there is a dense set of solutions of these 
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equations) and the integration over the p's, runs over the first Brillouin zone (with exclusion of p — which only 
contributes to the static part of the correlation functions). 
We will now focus on the density, which reads (for large L): 

(n x (t)) - (nx(oo)) = «««(«>)) - (»«(0))) / ^~ ipS ~ 6pt {f Pl - eiP(2 ~ l) f f*(a,j) ) ' {P * 0) (W) 

With the boundary conditions, we obtain for the stationary density: 

i) For the SEP model, where 03 + 05 = o^dp.o'B (i- c - |B| = C + -D) with o- 2 = og > and o\ + 04 = A = 0, we 
recover the results of Stinchcombe and SchutzaQ. 

ii) For the other cases on the 10— parametric manifold (Q), for which \B\ ^= C + D, i.e. 0-3 + 0-5 7^ er 2 + °6, and thus 
A = 0-1 + 0-4 7^ : 

, / x\ , . -V2T4V + t 2 (t 3 -(J 6 f] L - 2 ) t 4 - o- 6 ?7 L - 2 ((rii(oo)) - (f) 

(r/o- 2 +T 1 )(r 3 +o- 6 77 L - 2 )- ( j 2 o- 6 77 i - 1 r t 3 + a 6 ?] L - 2 

(n x ((x>)) = (( ni ((x>)) + (n L ((x)) - 2<j))r 1 x - 1 + (f), I < x < L, (15) 



where, ef> ee ^jg^ i ^ 2±^g^El^ ; n s (1 _ 2 ( Q + j))a 6 (a 3 + a 5 )^ r 2 ee 2a 4 - 5 + 
[2(o- 6 - 0-5) + a + 7] ; r 3 ee (1 - 2(/3 + (5))o" 2 - (0-3 + 05); and t 4 ee 2<7i - S + cf) 2(cr 2 - 03) + /3 + 5 
The dynamical part of the density reads, according to (p"4|), 

i) For the SEPjupdel, where B + C + D = 0-3 + 0-5 — 0-2 — o-g =0 and A = 0-1 + 0-4 = 0, we recover results obtained 
in references Rcf.pij. 

ii) For the other processes on the 10— parameter manifold (Q) (for which 03 + 0-5 ^ o^ + og with A ^ and o^og 7^ 0), 
we have 

L 



(n x (t)) - K(co)) = -e Bt J2 (§) 2 «»»(0)) - (%(oo)» (W2|£|*) + AJ x+J ,_i(2|E|t)) 

£ y ^ X—k 

+ e s *(l-A 2 )X:^ A ^ 1 f??) 2 (K(0))-K(oo)))4 +!/ _ 1 (2|S|t), (16) 



where A ee ,/a(l - 2(a + 7)) = 2 *?7 (a+7) . 

Similarly, we obtain the expression of the non-instantaneous two-point correlation functions (n x (t)n Xo (0)} for an 
initially translationally invariant system with initial density (n m (t = 0)) = p(0), m = 1, . . . , L. In so doing, we have 
to replace in (fil| (n.j(0)) with /O 2 (0)(l - 5 Xja;o ) + ^(O)^^. 



IV. THE QUENCHED DISORDERED SITUATION 



In the subsequent sections we will exploit the results of the previous section to compute the density and non- 
instantaneous correlation functions of systems obeying ([!]) and in presence of quenched disordered. More specifically, 
we assume that the reaction-rates defined in (||) obey a specific probability distribution. We assume a distribution, 
which has already been considered in the study the Glauber- Ising Model and other diluted chains (see e.gBtSEJ), and 
more recently, of the symmetric one-dimensional single-species diffusion-limited reaction A + <-*■ + An. The 
first distribution under consideration is the so-called "peaked distribution" which will describe a diluted chain: 

V(A,B,C,D) = tf k£fi 6 A ,A 1 &B,B 1 &c,Ci&D,iH + (1 - q &C )Sa,oSb,oS c ,o#d,0) where < q < 1 and A, B,C, D have 
been defined in (|^). According to this distribution, the chain is broken into discontinuous segments of (finite) length 
C, with C + 1 sites. On these segments the reaction-rates are respectively A%, Bi, C\, D\. Therefore, each realization 
has a probabilistic weight lo(C): 

u(C)=q &c \&J>0, (17) 
if the segments under consideration have periodic boundary conditions and a probabilistic weight u>(C): 
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lo(C) = J (1 " q&)2q ' * 1 " C ~ L ~ 2 fl8) 
1 (1 - g*)^' , if £ = £ - i 

if the segments under consideration are open segments. 

In Rcffl, such a distribution (18) has been considered for the symmetric A + «-> + ^4 reaction, for an initially 



homogeneous and uncorrelated system. The authors used the relationship of their model to the Heisenberg chain and 
its symmetries to compute explicitly the autocorrelation function (n r (t)n r (C 
For periodic segments, we also consider the algebraic disorder-distributior 

(n _ J , if C < b 

[ ' ~ I (7 - l)b^C^ , if C > b; with 7 > 2 Uyj 

In this section we want to extend Grynberg and Stinchcombe'si results and calculate explicitly the density and the 
non-instantaneous two-point correlation functions on the 10— parameter space manifold (|l]) on which the equations of 
motion of the correlation functions are closed. This is first performed for disordered systems of periodic segments of 
length C and then for disordered systems of open segments. 

To obtain the density of the disordered systems, which we denote (n m )(t), and of the non- instantaneous two-point 
correlation functions (n m (i)n;(0)), (where the symbol (. . .)(t) means the average over the time and the disorder), we 
need the density (n m (t))c and the non-instantaneous two-point correlation functions (n m (t)ni(0)} £ for a finite and 
ordered system of length C. In this section we consider the disorder-averaging of an observable, say O, as: 

We E^(W)) £| (20) 

where C = J £2 c ui(£)£, and (0(t))c denotes, respectively, the (timc-)mean- value of the segment of length L and 
the (time-)mean-value of an observable O on a lattice of size L. In (|20|), the averaging over the samples (i.e. the 
quenched- disorder) is carried out over segments of length C. It has to be stressed that this summation must be 
compatible wit h the na ture of the observable O. In fact, if one is interested in computing the density of particles at 
site m, namely (n m (t)}, one should carry the summation over all the segments which include the site m. In this case, 



(n m (t)) = — — — and similarly for the correlation functions. 

To study the mean value of an observable O, which time-average on a finite segment of length C is (0(t))c, we 
distinguish two regimes: 

To study the quantity (0(t)) (pp|), we will distinguish two long-time (i.e. \E\t ^> 1) regimes: the regime where 
u = L 2 /\E\t ~ 1, i.e. for typical diffusion-time and the regime where u = L 2 /\E\t <C 1, where time is much larger 
than the typical diffusion-time. In the latter the finite-size effects due to the disorder sampling procedure dramatically 
affect the dynamics. There is also a third regime, where |£7|t 3> 1 and u = L 2 /\E\t 3> 1, which is less interesting for 
our purpose, where we have (0(t)) ~ (O(t))oa. In the sequel we will focus on the regimes u ~ 1 and u -C 1 and omit 
the case u ^> 1. 

We are interested in the long-time behaviour of statistical quantities such as density or correlation functions (see 
section II. A) and (|l6|). It has been established (see e.g. Ref.Bi3) that the long-time behaviour of these quantities 
are governed by the small momenta p, ergo by large segments, with 1 <C C ^ L — > oo. Here the observables 
O considered are the density operator at site j, i.e. O = rij and the operator UjUk- For these operatorSj— as 
discussed, the long-time behaviour on a finite segment of length C follows from a small momentum expansionlij : 
(0(t)) - e ( B + 2 l B D' J2 P ex P ("l^l^ 2 (! + 0{l/C 2 ))) , according to the results of section II.A and © The long-time 
behaviour (\E\t 3> 1) in the presence of quenched-disorder, in the regimes u ~ 1 and u< 1 follows aa3 : 

,if_l 

^c <V^Tmt , if u « 1 



Let us note that the SEP model with free boundary conditions, i.e. a — (3 — ^f — 5 — 0, is special in the sense that 
in this case the quantization constraint ( [l3|) of the momenta p's is always ifor segments of arbitrary length £, and 
not only for the case C w L — > oo) fulfilled with p = , j = 0, 1, — IB. 
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V. QUENCHED DISORDER FOR PERIODIC SEGMENTS 



In this section, we compute the density and the non-instantaneous two-point correlation functions on the soluble 
10 -parameteric manifold for systems with periodic boundary conditions ^. 

For these periodic systems, we assume that quenched- disorder is characterized by two different distributions: the 
exponential distribution ([l7]) and the power-law distribution (19). 

1) We begin with the first distribution ([l7]) and introduce the following notations: £ = — s = j. 
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^ ; A=a^(^y +2 +2^e^(^y /+2 ; andS^d/3(/3-l)(^)^ + 12^(^)"^ 

l.i) To study the density and the correlation functions in this regime, we first expand the expressions obtained in 
section III. A at small momenta (large C) and then proceed to the disorder-averaging according to ( |2~i| ) within a saddle 
point expansion. 

In the regime where u < 1, we have for the density of a system with a non-uniform initial density profile : 
(n m (t = 0)) = (n m (t — oo))(l + 5 m ,i)- With r = m — 1, we have: 



/ 2^(n m (oo))e re *\ 11 



(n m (t)) - Moo)) « x T'' * p= '- (22) 

For the non- instantaneous correlation functions, we also consider 

a non uniform initial system: (n y (0)n Xo (0)) = (n v (oo)n Xo (0))(l + 5 Vt i)(l — S VtXa ) + (n Xa (0))S y ^ Xa . The expression 
of (n x (t)n Xo (0)) — (n r (oo)no(0)) with such initial state is obtained from r.h.s. of ( p2|) in replacing (n m (oo)) with 
fx*-] (Moo)) +^-*«((n Xo (0))) ~ (1 +4 ,i)K(oo))). 

l.ii) In the regime where it ~ 1, according to (J2l|), we have for the density 

s 2-0 



M*)) - Moo)) a ( | ) ' (^l^lt) 1 7 e"f <* «n ro (i)> - (rc m (oc)» (23) 



For the non-instantaneous correlation functions, according to (21), we replace in the r.h.s. of (|23]) ((n m (t)) — (n m (oo))) 
with ((n r (t)n (0)) — (n r (oo)no(0))) Expressions (n m (t)) — (n m (oo)) and (n r (i)n o (0)) — (n r (oo)) have been obtained 
previously (|4j) and (0) for the (infinite) ordered periodic chain 

We now consider the second (algebraic) distribution ( |l9| ) for the disorder average, with 7 > 2. 

2.i) For a system characterized by a non-uniform initial density: (n m (t = 0)) = (n m (oo))(l + S m> i), with r = m — 1 



,{B+2\E\)t 



Mt)) ~ Moo)) « ^^-*( 7 - i)^r 1 j (47r2|g|f)(7 _ 1)/2 (24) 

For systems with non-uniform initial states such that (n y (0)n Xo (0)) = (n y (oo)n Xo (0))(l+5 y ^i)(l — S ytXo ) + (n Xo (0))6 ytXo , 
we have to replace in the r.h.s. of ( p4] ) (n m (oo)) with fi x ~ 1 ((n Xo (oo)) + ^ 1 ~ :Eo ((n :Eo (0))) — (1 + S X(u i)(n Xo (oo))) . 
2.ii) In the regime where w ~ 1, we have for the density, according to (|2l|), 



MQ) ~ Moo)) a ^- 1 (4 7 r 2 |£;|t) 1 -* «n ro (t)) - (n m (oc)» (25) 



For the non-instantaneous correlation functions, according to (|2lj), we obtain: (n r (t)no(0)) — (n r (oo)) 
^- 1 (47r 2 |i?|i) 1 -* «n r (i)n o (0)) - (n r (oo)n (0))). 
In these expressions (n m (t)) — (n m (oo)) and (n r (t)no(0)) — (n r (oo)) have been obtained in section III. A. 



2 For segments (with periodic boundary conditions as well as for open segments) caution is required in order to preserve 
the parity of the processes when considering the quenched-disordered steady-states of the system. As an example of processes 
which preserve the parity of the number of particles, such asA + $^$ + A and 71 + ^4^0 + or the process A + «-» + A 
and A + A — > + 0, one has to take into account that for (finite, and infinite) segments with initially odd number of particles 
there wili afways remain one particle on each segment (see Ref.B). That is the reason for the steady-state (e.g. (n r (oo))) to 
be considered separatly. In order to avoid the difficulties due to residual particles, we can restrict to the even subspace, i.e. 
^ (1 + PJ o"™)|P(0)). Therefore, we can, e.g., only take into account segments with even number of particles and the disordered 
steady-states correspond to (. ..)(< = 00) = = ^2 C lu(C)C(. . -)c(t — 00) 
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VI. QUENCHED-DISORDER FOR OPEN SEGMENTS 



In this section (with the same notations as in sections V and VI), taking advantage of the results of the section IV, 
we compute density and correlation functions on the soluble 10 -parametric manifold for systems with open boundary 
conditions. 

We consider the effect of quenched-disorder for open segments. We assume that the system is initially uncorrelated, 
with density p(0). For technical convenience, we assume the exponential distribution (|l8|). In what follows, the case 
of (finite) segments with free-boundary conditions (without injection and evacuation of particles at the sites 1 and C) 



corresponds to the case where a = f3 = ~/ = 5 = and thus a = 
We introduce the quantities ^ = ^(&,/3) , z = 1,2,3, where 



,P = h + 



2cr 2 



&=ln|A|±ln 



^3 = In h| ± In 



,7 



fT4— &6 X _ CTj 

2cr G ' U 2cr 2 ' 



(26) 



where, as previously, we consider D/C = o§jo~2 > and assume A 



2tr 5 -(a+7) 
2<r 6 



> 0, for a 5 > and DC = a 6 a 2 > 0, 



without further restrictions, we can always assume that 2o§ > a + 7 > 0. In particular, this holds true for free- 
boundaries systems, where a — [3 — 7 = 5 = 0. 

As in the previous section, we distinguish two regimes: 

i) In the first one, L 2 /\E\t <C 1, the result is obtained via a saddle-point expansion. 

To study the density and the correlation functions in this regime, we first expand the expressions fllq) and its 
correlation functions for small momenta (large C) and then proceed to the disorder- averaging according to (pl|) within 

a saddle point expansion. In so doing, the saddle points sq = So(£j) , z = 1,2,3 are solutions of the saddle-point 

1/3 

equation: - &/3(s £) /l ~ 1 + ^ = , i = 1, 2, 3 For the case $ = 1, we have s (£-) = (j=gr) ■ 
The long-time behaviour of the density after averaging over the time and the disorder reads (for 



(n*(*)> - M°°)> * ^ 



exp 



(B + 2\E\)er 
2\E\ 



max 



{«} 4=1>2 ,3 



2tt 2 t 

4m 



(27) 



where v = v^) is an exponent which depends on the parameters . On the r.h.s. of (|27j), in the term 
maxr f ,\ {■■■}, we select the maxima (over the £■) of the function /(so(£0) = ££; s o(£D ~ Q ; C^ _ls o(Ci) /3 — 

l**Jl=l,2,3 

yy, (i = 1, 2, 3), where the So(CD are ^he solutions of the saddle-point equation. 



mo 



We see from the formula ( |27| ) that the effect of the open boundary appears through the term £| ( |26| ) and can affect 
the dynamics with respect to the case of segments with periodic boundary conditions. 
For the case where $ = 1 all the computations can be explicitly carried out as follows. 



We define £' according to {a — 



2/3 



{f'}>=i 2 3 ~ ££;} 2//3 ■ To compute the density, we then distinguish 



two cases (with a ^ ££'). First, if £' ^ (or f = and A = 1, as for the SEP model). 



(n x (t)} - (n x (oo)) « — 



27T / 4?Ht 

J { (a - tfr 



1/6 



exp 



3(a - ft') 



/\2/3 



-(47r 2 r) 



1/3 



5 + 2|£| 
2|£l 



(28) 



If £' = 0, A 7^ 1 and f3 = 1, for the situation considered here (systems initially uncorrelated with initial density 
p(0) ^ 0), we have 



(n x (i)) - (n x (oo)) 



a 1 / 3 



• exp 



-(4 7 r 2 a 2 r) 1 / 3 + 



+ 2|£| 

2|£| 



(29) 



ii) In the other one, where L 2 ~ at /3 = 1, with a 7^ according to (|2l]), we have 

(Mt)) - KH) 



- K(oo)) pb -= (tt 2 ^!*) 3 exp 



-l(^l«lt)i 



(30) 



In this expression, (n x (t)) — (n x (oo)) is given in (|l6|), for the case |2?| ^ C + D. 

Similar computations can be performed for (n x (t)n Xo ) — (n x (co)n Xo ) . For systems initially uncorrelated 
and homogeneous with initial density p(0) and for the case (3 — 1, the long-time behaviour (\E\t ^S> 
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1, u 



L 2 /\E\t <C 1) of the two-point correlation functions is the following: (n x (t)n Xo ) — (n x (oo)n x 



exp 



3(a-£t 



-(47r 2 r 



,1/3 



g+2|Sj 
2|B| 



£ 2 t J , where £' is the quantity defined above. We see that in absence of 



initial correlations, the two-point correlation function has, in this regime, the same stretched-exponential long-time 
behaviour as the density (p8|). The exponent v(C) depends on For the cases considered here — 1) its value is 



i (as for the SEP model with free boundary conditions, see below and Ref.H) or v 
In the other regime, the long-time behaviour (\E\t ^> l,u = L 2 /\E\t ~ 1) of the two-point correlation functions 
follows as in @):(n x {t)n Xo ) - (n x {oo)n Xo ) ra ^ (n 2 \E\t) ' exp [-f (^ 2 |£|i)5] ((ra^K^O)) - (n a (oc)n Xo (0))) . 

In the presence of initial correlations, the situation changes radically: one should take into account the initial 
correlation in the disorder-averaging. The latter give rise to results which would be very different from the ones 
obtained for the density. 



VII. ILLUSTRATION: TWO EXAMPLES 



As a first example we consider the symmetric exclusion model (SEP) A + <-> + A, with rates rj" = r° J = T > 0, 
on a diluted chain (periodic and with free boundaries) with distribution probability: V(T) = qdr^t + (1 — q)$r,o- 
With the results of the section III. A, III.B, V and VI, setting a = j3 = 1, in the expression of A and B (see section 
V), in the regime Ft 3> L 2 , we recover the result of Grynberg and Stinchcombea: 



(n x (t)n X0 ) - (n x (oo)) « p(0)(l - p(0)) J^(4tt 2 t) 1 / 6 exp 



--(4^/3 



(31) 



This result is valid for open segments as well as periodic segments. 

As second illustration, we consider the dynamics of the system described by the following processes: A + A «-* + 0, 
with rates rJJ = = E and A + <-> + A, with rates T?J = Tj? = V. „ 

This model is relevant to the description of dimer adsorption and desorption (see Ref.cl, and references therein). 

Contrary to the SEP model, here the effect of open boundaries is real. We assume that Ti > Si and that 
there is neither injection nor evacuation of particles at the boundary (free boundary conditions) of a segment (i.e. 

a = = j = S = O). 

This model preserves the parity of the number of particles, a property which has to be taken into account when 
computing the steady states (see footnote 2). 

Consider a diluted chain described by the distribution probability: V(T, E) = qSr r^s Si + (1 — <z)^r ofe o > (0 < 
?<!)• 

Thus, for segments with periodic boundary conditions, the disorder-distribution simply reads ui(£) = q . 
If we assume segments with free-boundary conditions, the disorder-distribution reads: 
(1 - q) 2 q c , if 1 < £ < L - 2 
(l-q)q c , if£ = £-l 

It is useful to introduce the notations: 4> = \ ; < r\ = < 1 ; A = 1+ 2S * > ; B = -2(r 1 +E 1 ) ; E = 

r x - E x ; £i = ; and = In A ; C3 = m? 7 I fi = er 4 = E x ; cr 2 = cr 6 = Ii - E x ; cr 3 = er 5 = Ti + Ei ; 07 = 
2Pi ; er 8 = cr 9 = 2Ei 

We first consider the case of segments with periodic boundary conditions. For a system with non-uniform initial 
conditions: (n m (0)) = (n m (oo))(l +6 m 1) = |(1 + <5 m 1), we obtain (r = m— 1), in the regime where £, (Tj — Ei)t 3> 1 
, andTS fiwyt_ 

KW) " KM) « KM) (j^j /6 exp(^-^((2 7 ra) 2 r) 1 /3_iM^_^ (32 ) 

Let us now pass to an open system, with free boundary conditions, for which the initial state is uncorrelated with 
density (n m (0)) = p(0), (with p(0) ^= 0), we have: 

In the regime where |Ti |t, \Ei\t > L 2 > 1, with £' = In A > 0, (a ^ 

, 1/6 o ( 3(«-glnA) 2 / 3 



((^V) exp(-2J3=^l(4^r)Vs_^g) , if A < e 2 ^ 



<n,(t))-<n,(oo)>«< U-^-rz /V < * ' ' — _ (33) 



.^l/3y^r exp r 3(( 27r &)2 T ) 



"1 — £ 



if A > e 2 "/« 



We note on the basis of this simple example, the effects of the boundary, conditions. The SEP model, where the 
long-time dynamics does not depend on the boundary conditions (see RefJa) of the segments is a special case. 



9 



VIII. SUMMARY AND CONCLUSION 



We have considered the single-species one-dimensional models characterized by rates lying on a 10— parametric 
manifold. Such reaction-diffusion processes are (in arbitrary dimensions) formally completely soluble because the 
equations of motion of correlation functions are closed. After having generalized a previous approach to the dynamics 
of (ordered) systems with open boundaries, with help of the Dynamical Matrix Ansatz, we investigated the dynamics 
of this class of processes in the presence of quenched disorder (consisting in a random distribution of the reaction 
of periodic and open segments). We have computed the density and the two-point non-instantaneous correlation 
functions for various initial states. The study of the dynamics in the presence of quenched disorder leads to various 
regimes. The most remarkable effect of the disorder appears when time is much larger than the typical diffusion 
time. In this case the finite-size effect radically changes the dynamics: for the SEP model, the non-instantancous 
correlation functions no longer decay algebraically but follow rather a stretched exponential. For the other models, 
in this regime, the quenched disorder also affects the dynamics but less dramatically: the decay of the density and 
two-point correlation functions is a combination of the ordered exponential (corrected by subdominant terms arising 
from the disorder average) and a power-law arising from the quenched disorder average. Similar effects of disorder 
have been reported in recent experiments!! We also pointed out the non-trivial effect of the boundary conditions on 
the disorder average. We have illustrated our results by considering two physically relevant examples. 

The support of Swiss National Fonds is gratefully acknowledged. 
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The solution of this transcendental equations at small p leads to p = ±-7^ + 0(C 2 ). Therefore, it is sufficient to assume 
p = = 0, 1, ... in the large C expansion of the terms contributing to the density and correlation functions. 
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